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Abstract 

We investigate the imphcation of unparticle physics on the Bu^d {t^, K)tt decays under the 
constraints of the B^^s — B^^s mixing. We found that not only the unparticle parameters that 
belong to the flavor changing neutral current (FCNC) processes but also scaling dimension du 
could be constrained by the B^^s — B^^s mixing phenomenology. Employing the minimum 
analysis to the -Bn,d ~^ i^^j K)tt decays with the constraints of B^^s mixing, we find that the 
puzzle of large branching ratio for B^ and the discrepancy between the standard model 

estimation and data for the direct CP asymmetry of B^ K~^it^ and B^ — > vr^vr" can be 
resolved well. However, the mixing induced CP asymmetry of Bd Kstt^ could not be well 
accommodated by the unparticle contributions. 
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I. INTRODUCTION 



Recently some incomprehensible phenomena at B factories have been explored, especially 
Bu4 (tt, K)7t decays. Firstly, the observations on the large branching ratio (BR) for 
Bd ^ 7r°7r° decay with the world average B{Bd 7r°7r°) = (1.31 ±0.21) x 10"^ and the direct 
CP asymmetry for Bd tx^tx^ with Acp{Bd tt+tt^) = 0.38 ± 0.07 jl| are inconsistent 
with the theoretical estimations of around 0.5 x 10~^ and 10 — 20%, respectively. Secondly, 
a disagreement in the CP asymmetries (CPAs) for Bd — * K^ti~ and B^ K^ti^ has 
been observed to be —0.097 ± 0.012 and 0.050 ± 0.025 [l|, respectively, while the naive 
estimation is Acp = Acp{B^ — > 7r°) — ^cp(i?rf —> K^ti^) ^0. Although many 



theoretical calculations based on QCDF PQCD [3| and SCET j4| have been tried to 
produce the consistencies with data in the framework of standard model (SM), however, the 

n 

results have not been conclusive yet [5|. For instance, the recent PQCD result for the ^cp 
is 0.08 ± 0.09, which is actually consistent with the data. However, the PQCD prediction 
AcpiB'^ K^TT^)pocn = — 0.0llo;o5 still has lAa difference from the current experimental 



data 



I PQCD 



6[. In addition, the difference between (sin 2/5)^^,^0 and (sin 2/5) in the mixing- 
induced CPA from the PQCD prediction is 0.065 ± 0.04, which shows about 2cr off the data 
—0.30 ±0.19. Hence, the inconsistencies between data and theoretical predictions provide a 
strong indication to investigate the new physics beyond SM. 

There introduced many extensions of the SM, and enormous studies have been done 
on searching some specific models beyond SM, e^. on supersymmetric model j?!, extra- 
dimension model 8|], left-right symmetric model |9| and flavor- changing Z' model 10|. Al- 
though new physics effects will be introduced, however, in phenomenological sense we just 
bring more particles and their related interactions to our system. Recently, Georgi proposed 



completely different stuff and suggested that an invisible sector, dictated by the sea' 



11 



e m- 



12|. 



variance and coupled weakly to the particles of the SM, may exist in our universe 
Unlike the concept of particles in the SM or its normal extensions where the particles own 
the definite mass, the scale invariant stuff cannot have a definite mass unless it is zero. 
Therefore, if the peculiar stuff exists, it should be made of unparticles ll|. Furthermore, in 
terms of the two-point function with the scale invariance, it is found that the unparticle with 

n 

the scaling dimension du behaves like a non- integral number du of invisible particles |11| . 
Based on Georgi's proposal, the phenomenology of unparticle physics has been extensively 
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For illustration, some examples such as t u + 14 



studied in Refs. jll, 

and e~^e~ — > fi~^fJ'~ have been introduced to display the unparticle properties. In addition, it 
is also suggested that the unparticle production in high energy colliders might be detected 



by searching for the missing energy and momentum distributions 



12 



13|. Neverthe^ 



we have to point out that flavor factories with high luminosities, such as SuperKEKB 



n 



ess, 

3, 



SuperB [18| and LHCb [19| etc, should also provide good environments to search for the 
unparticle effects in indirect way. 

Besides the weird property of non-integral number of unparticles, the most astonished 
effect is that an unparticle could carry apeculiar CP conserving phase associated with its 



propagator in the time-like region 12, ll3||. It has been pointed out that the unparticle 



ay s a role like a strong phase and has an important impact on direct CP violation 



14( 1 . In this paper, we will make detailed analysis to examine whether the puzzles in 



phase p 
(CPV) 

Bu,d ^ (tt, 1^)7^ decays with the B^^s ~ B^^s mixing constraints could be resolved when the 
invisible unparticle stuff is introduced to the SM. 

In order to study the flavor physics associated with scale invariant stuff, we follow the 
scheme proposed in Ref. HI]. For the system with the scale invariance, there exist so- 
called Banks-Zaks {BZ) flelds that have a nontrivial infrared flxed point at a very high 
energy scale 20|]. Subsequently, with the dimensional transmutation at the Ku scale, the 
BZ operators composed of BZ flelds will match onto unparticle operators. We consider 
only vector unparticle operator in the following analysis. Then, the effective interactions for 
unparticle stuff and the particles of the SM are adopted to be 



a; 



7.(1 - iMOl, + -^9 7.(1 + , 

A7 



(1) 



where Cl'^j^ are effective coefficient functions and denotes the spin-1 unparticle operator 
with scaling dimension dy and is assumed to be hermitian and transverse d^Oj^ 



0. Since so 



far the theory for BZ flelds and their interactions with SM particles is uncertain, here 



L,R 



are regarded as free parameters. With scale invariance, the propagator of vector unparticle 
can be obtained by [12, ll3|] 



rf^xe'^-(0|T(O^(x)O^(0))|0) 



fiiy 



p2 



-i<f>U 



(2) 
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with = {du — 2)-7r and 



A 



167r5/2 T{du + 1/2) 



(3) 



{27^y^^^ T{du - l)T{2du) ' 

where (pu could be regarded as a CP conserving phase 1^, [l^. We note that when 
unparticle stuff is reahzed in the frameworlc of conformal field theories, the propagators 
for vector and tensor unparticles should be modified 2l|. Although conformal invariance 
typically implies scale invariance, however, in principle it is not necessary. Unparticle stuff 



221. Hence, in this 



with scalar invariance in 2D spacetime has been investigated in Ref. 
work, we still concentrate on the stuff built out of only scale invariance. For simplicity, we 
set the unknown scale factor Au to be 1 TeV throughout the analysis. 



II. CONSTRAINTS OF Bd,s - MIXING 

Since the tree level FCNC processes are allowed in the unparticle physics, it is expected 
that the Bd^s — Bd^s mixings offer strong constraints on the unparticle parameters of Cf', 
and C£*, C|j*. Moreover, it will be shown that the scaling dimension du also can be con- 
strained by Bd^s — Bd.s mixings. 

First of all, we separate the matrix element of /S.B = 2 transition for the 5° — B^ mixing 
{q = d,s), which is denoted by MI2, into the SM and unparticle contribution as follows. 

Mf, = Mli^' + Mff ^ = \Ml^''''\e<"' + |Mff ^|e<" . (4) 

where the (p^^ and (p^^ represent the phases of mixing amplitudes. For the second term, 
we use the superscript 'NP' in order to represent general new physics (NP) contribution. 
Later on, we regard this Ml^^ as the unparticle mixing amplitude. 

As is well known, the magnitude of total mixing amplitude |Mj^2l is given by the B^ — B^ 
oscillating frequency as follows 

AMg = 2\Mf2\. (5) 

And the mixing phase 0g = arg can be obtained from the mixing induced CP asymmetry 
of 6 — > CCS processes. We summarize current experimental data in Table [B 



4 



TABLE I: Experimental values for the B^^s — B^^s mixings. Even though DO collaboration recently 



have measured 



23], the data is not used in our analysis because of huge error of it. 



observables 


values 


note 




( 0.507 ± 0.004 ) ps-i 


HFAG []J 


AM, 


( 17.77 ±0.12 ) ps-i 


CDF [24] 


4>d 


43° ± 2° 


HFAG [IJ 



The SM mixing amplitude reads 



M- 



q,SM 
12 



(6) 



where fj^ = 0.552 is short distance QCD correction term 25|, and Sn(xt) = 2.35 ± 0.06 is 
an Inami-Lim function for the t-quark exchange in the loop diagram 26|. The quantities of 
fsg and Bb^ are non-perturbative parameters which can be obtained from the lattice calcu- 
lations. We follow the procedure given in Ref. j^^] for dealing with these non-perturbative 
parameters. The procedure mainly employs the result of lattice calculations in two different 
ways. The one is to use the result of 



the results of JLQCD and HPQCD 



LQCD collaboration 



and the other is to combine 



291 ] collaborations. We note that one can obtain the 



SM mixing phases from Eq.([6]) as follows. 

2P, 



lSM 



iSM 



(7) 



where is an angle of CKM unitarity triangle, A and t] are from the Wolfenstein parametriza- 



tion 



301 ]. We use the result of UTfit [3l| from the tree level processes for the (3,Rt and fj, 
where Rt = a/(1 — p)^ + f]'^ and (p, f/) is the apex of the CKM unitary triangle. For the ]V^b], 
we adopt the result of global fit to moment of inclusive distributions in i? — > XJui, which 



is performed in the framework of heavy quark expansions with kinetic scheme [32|. After 
putting all the SM input parameters into Eq. the SM mixing amplitudes are obtained. 
And using the experimental data shown in the Table [H the NP mixing amplitudes for the 
Bd — Bd mixing could be gained through the Eq. (jl]) . The numerical values are summarized 
in Table M 

As for the unparticle contribution to the mixing amplitude, we begin with the effective 
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TABLE II: Numerical values for the SM B^^s — B^^s mixing amplitudes. The values of the NP 
Bd — Bd mixing amplitude are obtained from the experimental data and the Eq. (j3]). The case (a) 



denotes JLQCD, while the case (b) denotes (HP+JL)QCD. Ml. 



NP 



12 



Mf^ should be considered 



{q = s,d). 



SM parameters 


values 




NP parameters 


values 


2\M'^f'''\ 


' "-'-0.26 

0.97 ±0.29 ps-i 


(a) 
(b) 


2|Mff^| 


0.25 ±0.26 ps~i (a) 
0.46 ± 0.29 ps"^ (b) 




45.2° ±5.7° 






-130° ± 180° (a) 
-132° ± 12° (b) 




16.4 ±2.8 ps-i 
23.8 ± 5.9 ps-i 


(a) 
(b) 






-2.3° ±0.2° 







Hamiltonian for AB = 2 processes in unparticle sector such as 



q,U 



4 U^, 



du-1 



A. 



U 



X 



± 



p2 2 sin duT^ 

-11, (Cf (1 - 75) ± ± 75)) b qr (Cf (1 - 75) ± Ct{l ± 75)) h 

^-q i> (Cf (1 - 75) ± Ct{\ ± 75)) h q i> (Cf (1 - 75) ± Ct{l + 75)) h 



The factor 2 is from the fact that there are s and t— channel which give same result, and 
the factor 1/4 is due to the Wick contraction factor 33|]. From this effective Hamiltonian, 
the transition matrix elements can be shown as 



M' 



q,U 
12 



(A^) 



2 \<iu-^ 



(9) 



where the is defined by 



q,U 



qb\2 



5 rn\ 
12 p2 



± ctct 



5 
3 



6 p2 



(10) 



with = Tn?B^. A^(p^) is given in Eq. Q. In order to get the constraints on the unparticle 
parameters from the NP parameter regions shown in the Table HTl we first consider the phase 
of Mf^ . It depends on the scaling dimension du through the e"**^" term and the sign of 
sm^duTi) in the Ayi^p"^). The plot of arg(Mf2") versus du is shown in Fig[TJ 
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c3 -60° - 
-120° - 
-180° - 



1.0 1.5 2.0 2.5 3.0 

FIG. 1: Plot of arg(Mf^") versus du within the range of (-180°, 180°). 

It is interesting that axg{M'l^) takes only positive value as displayed in Fig. [H Therefore, 
as we can see from Table HTl (HP+JL)QCD case can not give a solution of du while all value 
of arg(M;^") is possible in JLQCD case. So, we take only the JLQCD case. Although we 
can not determine du in the JLQCD case presently, if the SM prediction for the Bd — Bd 
mixing amplitude becomes more precise in future, it can give strong constraint on du from 
the plot of arg(M{^2^) versus du- Here, we assume du = 1.5 for the remaining analysis. 
Then, the magnitude of unparticle transition matrix element M^^ can be obtained through 
Eq. (jlj) and the values in Table [III for the JLQCD case with the experimental data as 
2\M^^\ = 7.6 ± 6.6 ps~^. Using this value and the value for 2\Mf^\ in Table HTl we can see 
that the mixing parameter should be strongly suppressed as follows. 

\aZ\ = (1-1 ± 1-3) X 10-« , la^^J = (2.6 ± 2.3) x 10"^ . (11) 
And, the Eq. leads to 

\Ct-C^^\ = 2^^. (12) 

Therefore, the parameters Cf^ and (g = d, s) are strongly correlated under the Bd^s — Bd,s 
mixing. In the next section, these strong constraints on the unparticle parameters will be 
used for fitting the parameters in Bu,d — > (tt, K)n decays 
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III. Bu4 (vr, K)tt decays AND THE UNPARTICLE CONTRIBUTIONS 

I, the effective Hamiltonian for b —* qq'q' decays is 



According to the Eqs. ([T|) and 
obtained by 



u 



Cu{q') (Ct{qh)v-A + Ct{qh)v+A) (cp' {q'q')v-A + Cp' {q'q' 



IV+A , 



(13) 



where q = (d, s), q' = {u, d, s, c), {f'f)v±A = f'l^,i^ ± 75)/ and 



Cu{q' 



(14) 



Based on this effective Hamihonian, we study the unparticle contributions to the decay 
amphtudes for B^^d {t^i K)ti. It is known that the most uncertain theoretical calculations 
for two-body exclusive decays are the QCD hadronic transition matrix elements. To deal 
with the hadronic matrix elements, we adopt recent perturbative QCD (PQCD) calculations 
for the SM amplitudes and naive factorization (NF) approach for the amplitudes of unparticle 
contributions. The SM amplitudes for Bu^^ ^(^! -^)^ decays can be parameterized in the 



context of quark diagram approach (QDA) 34] as follows 



(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 

where T^'^ and C^'^ denote the tree color-allowed and -suppressed amplitudes for B^^ 
'K{K)7r, respectively, while P^'^P^ly) is gluonic (electroweak) penguin amplitude. All CP- 
conserving phases are included in these parameters. The phase 7(/5) is the CP violating 
phase in the SM and from Vyj, ( Vtd ) ■ T able UTTl shows the recent PQCD result for the values 



V2A''^\B+ - 


— i> tt" 


^7r°) -- 


= -Te'^ 


- Ce'"' - PEwe-'^, 


A''\Bd - 






= -Te*^ 


- Pe-'f^, 


V2A''\Bd 




TT ) - 


= -Ce'^ 


+ Pe-'f - PEwe-'^ 




^ K\+) -- 


= P\ 




A''\Bd - 




IT ) : 


= -P'- 




V2A''^\B+ - 


K~ 


*-7rO) -- 


= -P'- 


rpl i-i r^lil p' 


V2A''\Bd - 


K 


OyrO) - 


- P' n'p'^'y p' 



of each topological parameters 



35] 



For deriving the unparticle contributions, the definitions for relevant decay constants and 
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TABLE III: Recent PQCD predictions for the topological parameters of B^^d (t^, K)tt in unit 
of 10~^ GeV . The phases are indicating strong phases of the parameters in radian unit. The 
predictions include NLO calculation. 



Topology 


Abs 


Arg 


Topology 


Abs 


Arg 


P' 
T' 
P' 

a 


1 7+0.9 

-■^•'-0.6 


Z.J_Q 2 

0.1 ±0.0 
-1.3 ±0.1 
-3.0 ±0.0 


T 
P 
C 
Pew 


23.2^1° 
^•"-0.8 

4 3+2-1 

^- '-0.1 


0.0 ±0.0 

-1.1 ±0.0 
-0.1 ±0.0 



form factors are given by 



{P{p)\q75u\0) = ifpml^ 



{P{p)\qi,h\B{pB)) 



m 



(Pb+p),, T% 



m 



(22) 



with P = (vr, K), q = Pb — P and mp = m\)/{mq ± m^). Here, due to nip ^ m^, we have 
neglected the m?p effects m. B ^ P transition matrix element. Subsequently, by considering 
various flavor diagrams in which the typical diagrams mediated by unparticle are illustrated 
in Fig. [2], the unparticle amplitudes for Bu^ (vr, K)tt decays within NF approach are 

s, d ^ ^ q 



U_ 




s, d 



FIG. 2: Typical flavor diagrams for B (vr, K)^ decays mediated by unparticle with q=u and d, 
where qi^2 are the momenta of unparticle. 



obtained to be 



Cu{ql)UrnlFt{ml)a'i::-' 



Cuiqi)fKmlFt{mj,)a'if^ 



(23) 
(24) 
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TABLE IV: The definition of coefficients for tlie unparticle amplitudes in -B^.d (tT) K)tt decays 
witliin NF approach. 



decay mode 



vr ' vr 



M 

-*dcc 



Nr \ \ L L R R J ' 1 \ L R R ^ 



T) 



V2Nc 



[Cf{C 



^ -ci^)-ci^{c^^-ci^)) 



(gj ) ( (jdh I (jdb\ / (juu (jdd (juu i (jdd\ 

y^Cuiqj) ^ ^ R ) \ L L R 



1 / / rjdbrjdd rjdbrjdd\ i '2f'^ ( fjdbrjdd rjdbrjdd'\ 

y^JSf^. \ \ L L R R J ' 2 \ L R R L J 

I (jd,b\ ( (juu (jdd (juu i (jdd\ 

R J \ L L R, R, J 



2T I'-^L 



J_ 



Irjsbrjdd rjsh rjdd\ i '2f^ (^sbrjdd ^sbr~<dd\ 

\ L L R R J 1 \ L R R L J 



Ifjsb^uu ^sb ^uu\ _j_ ^rpK /rjsbrjuu rjsb rjuu\ 

\ L L R, R J 1 \ L R R L J 



l^sb^juu ^sbf~<uu\ _j_ l^sb^uu ^sb^uu\ 

^y2Nc \ ^ L L R R J 1 \ L R R L J 

^u{<l2) fir ^Q^i'i^'i) I (~isb I r<sb\ ( r^uu /~idd /~iuu i ridd\ 



V2Cu{q^,)fKFB-(ml) 



(Cf + ) {Cf - cf - + cf) 



where the coefficients d^^^s are defined in Table HVl A^'c = 3 is the number of colors, 



mz 



and the chiral enhanced factor rj^ 2) 2) ^"^^ defined by 



m: 



mz 



K 



m 



K 



(25) 



We note that since qi in Fig. [2](a) involves different mesons, the estimation of q\ should have 
ambiguity. To understand the typical value of we write the qi = Pb ~ k2 — with /c2,3 
being the momenta of valence quarks inside the light mesons. In terms of momentum fraction 
of valence quark and light-cone coordinates and by neglecting the transverse momentum, 
one can get /c2 = {^,'^3X2/ V^,^±) and k^, = (m^xs/v^, 0, 0_l). As a result, we have 
qf = — 3^2) (1 — According to the behavior of leading twist wave function of light 

meson, <|)t^^2 oc x(l — x) which is calculated by QCD sum rules 37j, it is known that the 
maxima of ^2,3 occur at X2 = X3 ~ 1/2. Therefore, for numerical estimations, the momentum 
transfer could be roughly taken as qf ~ with ms = 5.28 GeV. 
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We discard irrelevant factor i in the NF and match the sign with the QDA parametriza- 
tion; then, the total amplitude is 

A(B ^ /) = KfA^'^iB ^ /) + A^{B ^ /) . (26) 

Here, the kj is the ratio of phase space factor coming from the difference of notation of 
decay amplitude between NF and PQCD group and is defined by 



-/-t/('^Vf5^)="5xlO-. (27) 



1287r J ' VSvrm^ 



From Table HVt we see that besides and A^, the introduced new free parameters are 

/^dfe /^db /^sb /^sb 

/^uu /^uu /^dd /^dd /'OQ^ 
1 1 1 ^R 1 \'^°) 

which denote the couplings of unparticle to SM particles. First four parameters are strongly 
correlated by Bd^s — Bd,s mixing phenomena as shown in Eg. (1121) . If we regard all these 
parameters to be real number and set Ku = 1 TeV and c/^ = 1.5 as we do in the Bd,s — Bd,s 
analysis, 8 free parameters are involved for Bu^d (tt, K)ti in the unparticle physics. 

In order to fit to the data for the Bu^d (^r, K)7r decays with these 8 parameters, 
we perform the minimum analysis. According to current experimental observations, 
the amount of available data for Bu^ (ti", K)tt decays is 17 as their world averages are 
displayed in Table |V] Besides the BRs, the important quantities to display the new physics 
effects are the direct and mixing induced CPAs, where they could be briefly defined through 

^^ = TTp^' ^^ = rTp^' ^^^^ 

with A/ = e^'^'^Aj/Af, respectively. For direct CPA, / could be any possible final states; 
however, for mixing induced CPA, / could only be CP eigenstates. Since the PQCD pre- 
diction for the SM decay amplitudes has sizable error as shown in Table llllt it should be 
considered for the minimum analysis. Therefore, we define the to be 

Ci and ti denote the experimental data for i'th observable and its theoretical prediction within 
unparticle contribution, respectively, af^^ is the experimental error of z'th observable, while 
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TABLE V: The experimental data for B^^a (tt, K)tt decays and comparison between the 
experimental data and the theoretical predictions with and without unparticle contribution. The 
BRs are order of 10^^. The data is updated by September 2007. 'w/o' means 'without unparticle 
contribution', contributions of each observables are shown. 



observables 


data 


theory (w/o) 


theory 


(w/o) 

A, \ / / 




6(/i°7r+) 


23.1 ± 1.0 


23.5 ± 12 


23.1 ± 11 


0.001 


0.0 




12.9 ± 0.6 


13.0 ± 6.2 


12.7 ± 6.0 


0.001 


0.001 


B(K+TT-) 

\ / 


19.4 ±0.6 


19.7 ± 10 


20.3 ± 10 


0.001 


0.007 


B(K'^TT^) 


9.9 ±0.6 


8.8 ±4.9 


9.5 ±5.1 


0.046 


0.006 




0.009 ± 0.025 


0.0 ± 0.0 


0.0 ± 0.0 


0.13 


0.13 




0.050 ±0.025 


-0.017 ±0.068 


0.074 ± 0.068 


0.84 


0.11 




-0.097 ±0.012 


-0.099 ± 0.073 


-0.11 ±0.075 


0.001 


0.03 




-0.14 ±0.11 


-0.065 ± 0.040 


-0.058 ± 0.036 


0.41 


0.50 




0.38 ±0.19 


0.74 ±0.08 


0.74 ±0.07 


3.1 


3.2 


H(7r+7rO) 


f •OJ-0.40 


4.03 ±2.53 


4.05 ±2.53 


0.37 


0.36 


B{tt+tt~) 


5.16 ±0.22 


6.80 ± 4.43 


7.11 ±4.43 


0.14 


0.19 


Bin^TT^) 


1.31 ±0.21 


0.23 ±0.13 


1.33 ±0.30 


19 


0.002 


Acp{7r+7T^) 


0.06 ±0.05 


0.00 ±0.01 


0.055 ±0.018 


1.6 


0.01 


^C'p(7r+7r^) 


0.38 ±0.07 


0.17 ±0.10 


0.38 ±0.14 


2.8 


0.0 




0.48t°i? 


0.64 ±0.23 


0.53 ±0.13 


0.17 


0.024 


c 


-0.61 ±0.08 


-0.55 ±0.44 


-0.55 ±0.42 


0.021 


0.023 



cr* is theoretical error propagated from the errors of topological parameters obtained from 
PQCD. Since there are 8 free parameters involved in our analysis, the degree of freedom 



(d.o.f) for the fitting is 



^17 — 8) = 9. As for the angle 7, we use the values of 7 = (63 



+l5\o 
-12) 



from the PDG 2006 36[. Consequently, by imposing the mixing constraints of C^(^) and 



L{R) 

Cf^f^j^^^ displayed in Eq. ([T^ . we find the optimized values of unparticle parameters as follows: 



Cf = 3.3 X 10-^ 



4.6 X 10"^ 



Cf = 7.6 X 10"^ C]^ 



R 



11.2 X lo-^ 

11.2 (31) 



5.0, C]^" = 12.0 , = 4.2, 

with = 4.6, compared to = 28.8 without the unparticle contributions. From above 
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TABLE VI: Numerical values of unparticle amplitudes for each decay mode. 'Abs' represents 
the magnitude of the amplitude in unit of lO^^GeV. 'Arg' represents the CP conserving phase of 
unparticle in radian unit. 



decay mode 


Abs 


Arg 


decay mode 


Abs 


Arg 




18.9 


-1.6 


vr+vrO 


0.6 


1.6 


K+TT^ 


10.8 


1.6 


TT~^TT~ 


3.5 


1.6 


K+TT- 


2.2 


-1.6 




9.9 


1.6 




2.9 


1.6 









results, we see clearly Bd^s — Bs^d mixings give strict constraints on C^5^) and Cf(^) that lead 
to FCNCs at tree level. We compare the experimental data to the theoretical predictions with 
and without unparticle contributions in Table IVl in detail. And also, the numerical values of 
unparticle contributions are given in Table |VT] for comparing with the SM contribution given 
in Table IIIII Strikingly, we can see that the large experimental data of BR for Bd n^ir^ 
can be quite well accommodated with unparticle contributions. Moreover, some anomalous 
observables of direct CPA of B^ f^+vr^ and Bd tt+tt" can be explained. However, 
the puzzle of mixing induced CPA of Bd Kst^^ could not be resolved well by unparticle 
contributions. As many authors argued, sizable non-SM weak phase is required in order to 
fit to the data of Sk^ttO 
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Since the unparticle contributions do not carry any extra weak 



phase, it turns out to be very hard to fit to the data. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have studied the effects of unparticle on Bu^ decays with Bd,s — Bd,s 
mixing constraints. For simplicity, we concentrate on vector unparticle and set the scale 
of unparticle operator to be 1 TeV. With the lattice QCD results of JLQCD and 
(HP+JL)QCD on the non-perturbative quantity of /^^ \J Bb^, we start to calculate the 
SM predictions for Mff^^. Accordingly, from the current data we extract the available 
space for the new physics effects in model independent way. When the unparticle effects are 
included to AB = 2 processes, we find that the current experimental data for AMd, AMs 
and (pd could give strict constraints on unparticle parameters as well as scaling dimension 
di(. The (HP-(-JL)QCD case could not give a solution for du, while all value of du is possible 

13 



in JLQCD case. However, we see that more accurate SM prediction for the — mixing 
amphtude in future would give strong constraint on du- In order to understand whether 
the unparticle effects could satisfy all measurements in exclusive B^^ — > (tt, K)tt decays, we 
utilize the minimum analysis to search for the solutions of free parameters. Interestingly, 
after we fix = 1.5 for the specific unparticle scaling dimension, we find that the unsolved 
problem of large BR for Bd tt^tt^ could be explained excellently in the framework of 
unparticle physics. Moreover the discrepancy between the standard model estimation and 
data for the direct CPA of B'^ K'^ir^ and B^ tt+tt^ could be reconciled very well. 
However, the puzzle of the mixing induced CPA of Bd Ksn^ could not be resolved well 
in unparticle physics. 

Acknowledgments 

This work of C.H.C. is supported by the National Science Council of R.O.C. under 
Grant No. NSC-95-2112-M-006-013-MY2. C.H.C. would like to thank Dr. Shao-Long 
Chen for useful discussions. The work of C.S.K. was supported in part by CHEP-SRC 
and in part by the KRF Grant funded by the Korean Government (MOEHRD) No. 
KRF-2005-070-C00030. The work of Y.W.Y. was supported by the KRF Grant funded by 
the Korean Government (MOEHRD) No. KRF-2005-070-C00030. Y.W.Y thank Jon Parry 
for his comment on mixing constraints. 



[1] E. Barberio et al. [Heavy Flavor Averaging Group (HFAG) Collaboration]. l arXiv:0704. 35751 



[hep-ex], online update at http://www.slac.stanford.edu/xorg/hfag 



[2] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B606, 245 
(2001) | arXiv:hep-ph/0104110| ; M. Beneke and M. Neubert, Nucl. Phys. B675, 333 



(2003) [arXiv:hep-ph/0308039| ; X. q. Li and Y. d. Yang, Phys. Rev. D72, 074007 (2005) 
|arXiv:hep-ph/0508079] . 

[3] H. n. Li and H. L. Yu, Phys. Rev. D53, 2480 (1996) f arXiv:hep-ph/9411308| ; Y. Y. Keum, H. 
n. Li and A. I. Sanda, Phys. Lett. B504, 6 (2001) |arXiv:hep-ph/0004004| ; Phys. Rev. D63, 



14 



054008 (2001) [arXiv:hep-ph/0004173] ; H. n. Li, S. Mishima and A. I. Sanda, Phys. Rev. D72, 
114005 (2005) | arXiv:h ep-ph/0508041|. 

[4] C. W. Bauer, D. Pirjol, I. Z. Rothstein and I. W. Stewart, Phys. Rev. D70, 054015 (2004) 
|arXiv:hep-ph/0401188] ; C. W. Bauer, I. Z. Rothstein and I. W. Stewart, Phys. Rev. D74, 
034010 (2006) [arXiv:hep-ph/0510241| ; A. R. Wilhamson and J. Zupan, Phys. Rev. D74, 
014003 (2006) [Erratum-ibid. D74, 03901 (2006)] iarXiv:hep-ph/0601214| . 

[5] A. Jain, I. Z. Rothstein and I. W. Stewart. [arXiv:0706.339"9l 

[6] H. n. Li, S. Mishima and A. I. Sanda, Phys. Rev. D 72, 114005 (2005) [arXiv:hep-ph/0508041] . 

[7] C. H. Chen and C. Q. Ceng, Phys. Rev. D66, 014007 (2002) [arXiv:hep-ph/0205306| ]; S. 
KhaUl, Phys. Rev. D72, 035007 (2005) [arXiv:hep-ph/0505151| ; R. Arnowitt, B. Dutta, B. 
Hu and S. Oh, Phys. Lett. B633, 748 (2006) [arXiv:hep-ph/0509233| ; Phys. Lett. B641, 
305 (2006) |arXiv:hep-ph/0606130| ; C. H. Chen and C. Q. Ceng, JHEP 10, 053 (2006) 
|arXiv:hep-ph/0608166| ; K. Cheung, S. K. Kang, C. S. Kim and J. Lee, Phys. Lett. B 652, 
319 (2007) f arXiv:hep-ph/070205 0l. 

[8] K. Agashe, G. Perez and A. Soni, Phys. Rev. D 71, 016002 (2005) [arXiv:hep-ph/0408134] ; 
S. Chang, C. S. Kim and J. Song, JHEP 0702, 087 (2007) [arXiv:hep-ph/0607313] . 

[9] R.N. Mohapatra and G. Senjanovic, Phys. Lett. B79, 283 (1978). 
[10] P. Langacker and M. Pliimacher, Phys. Rev. D62, 013006 (2000); V. Barger, C. W. Chiang, P. 



Langacker and H. S. Lee, Phys. Lett. B598, 218 (2004) |arXiv:hep-ph/0406126| ; D.A. Demir, 
G.L. Kane, and T.T. Wang, Phys. Rev. D72, 015012 (2005) |arXiv:hep-ph/ 0503290i ; C. H. 
Chen and H. Hatanaka, Phys. Rev. D73, 075003 (2006) |arXiv:hep-ph/0602140| ; S. Baek, 
J. H. Jeon and C. S. Kim, Phys. Lett. B 641, 183 (2006) |arXiv:hep-ph/0607113] . 



[11] H. Georgi, Phys. Rev. Lett. 98, 221601, (2007) [arXiv:hep-ph/0703260| . 
[12] H. Georgi, Phys. Lett. B650, 275 (2007) |arXiv:0704.2457] [heo-ph]]. 

[13] K. Cheung, W. Y. Keung and T. C. Yuan, Phys. Rev. Lett. 99, 051803 (2007) [ arXiv:0704.2588] 
[hep-ph]] ; larXiv:0706.3T55l [hep-ph] . 

[14] C. H. Chen and C. Q. Ceng, Phys. Rev. D76, 115003 (2007) [arXiv: 0705. 0689) [hep-ph]]: Phys. 
Rev. D76, 036007 (2007) |arXiv:0706.0850l [hep-ph]]: larXiv:0709.0235l [hep-ph]. 

[15] K. Cheung, W. Y. Keung and T. C. Yuan, l arXiv:0710.2230] [hep-ph]; larXiv:0711.336T] [hep- 
ph]; M. Luo and G. Zhu, Phys. Lett. B 659, 341 (2008) [arXiv:0704.3532] [hep-ph]]: M. x. Luo, 
W. Wu and G. h. Zhu, Phys. Lett. B 659, 349 (2008) [arXiv:0708.067T] [hep-ph]]; G. J. Ding 



15 



and M. L. Yan, Phys. Rev. D 76, 075005 (2007) [arXiv:0705.0794] [heo-ph]]: [aFXiv:0706.0325] 
[hep-ph] : la?Xiv:0709.3435l [hep-ph] : Y. Liao, Phys. Rev. D 76, 056006 (2007) [ arXiv:0705.0837] 
[hep-ph]]; larXiv:0708.3327] [hep-ph]; Y. Liao and J. Y. Liu, Phys. Rev. Lett. 99, 191804 
(2007) [ arXiv: 0706. 1284] [hep-ph]]; T. M. Ahev, A. S. Corneh and N. Gaur, Phys. Lett. B 
657, 77 (2007) [ arXiv: 0705. 13261 [hep-ph]]: JHEP 07, 072 (2007) [arXiv:0705.4542] [hep-ph]]: 
T. M. Ahev and M. Savci. [ajXrv:0710.1505l [hep-ph]: S. Catterall and F. Sannino, Phys. Rev. 
D76, 034504 (2007) [arXiv: 0705. 1664) [hep-lat]]; X. Q. Li and Z. T. Wei, Phys. Lett. B651, 
380 (2007) |arXiv:0705.182n [hep-ph]]; Xue-Qian Li, Yong Liu, Zheng-Tao Wei and Liang 
Tang, larXiv:0707.2285] [hep-ph]; C. D. Lu, W. Wang and Y. M. Wang, Phys. Rev. D 76, 
077701 (2007) [arXiv: 0705. 2909) [hep-ph]]: M. A. Stephanov, Phys. Rev. D 76, 035008 (2007) 
|arXiv:0705.3049] [hep-ph]] : P. J. Fox, A. Rajaraman and Y. Shirman, Phys. Rev. D 76, 075004 
(2007) |arXiv:0705.3092] [hep-ph]]: N. Greiner, Phys. Lett. B 653, 75 (2007) [arXiv:0705.3518l 
[hep-ph]]; H. Davoudiasl, Phys. Rev. Lett. 99, 141301 (2007) [arXiv: 0705. 3636] [hep-ph]]; 
D. Choudhury, D. K. Ghosh and Mamta, Phys. Lett. B 658, 148 (2008) [arXiv:0705.3637] 
[hep-ph]]; D. Choudhury and D. K. Ghosh, larXiv:0707."2074] [hep-ph]; G. Bhattacharyya, 
D. Choudhury and D. K. Ghosh, Phys. Lett. B 655, 261 (2007) |arXiv:0708.2835l [hep- 
ph]]; S. L. Chen and X. G. He, Phys. Rev. D 76, 091702 (2007) [arXiv: 0705. 3946] [hep-ph]]; 
S. L. Chen, X. G. He and H. C. Tsai, JHEP 0711, 010 (2007) [arXiv:0707.0187] [hep-ph]]; 
S. L. Chen, X. G. He, X. P. Hu and Y. Liao, larXiv:0710.5T29] [hep-ph]; P. Mathews and 
V. Ravindran, Phys. Lett. B 657, 198 (2007) |arXiv:0705.4599] [hep-ph]]; M. C. Kumar, 
P. Mathews, V. Ravindran and A. Tripathi. [a"rXiv:0709.2478] [hep-ph]; S. Zhou, Phys. Lett. 
B 659, 336 (2008) jarXiv: 0706. 0302] [hep-ph]]; R. Foadi, M. T. Frandsen, T. A. Ryttov and 
F. Sannino, Phys. Rev. D 76, 055005 (2007) [arXiv: 0706 . 16961 [hep-ph] ] : M. Bander, J. L. Feng, 
A. Rajaraman and Y. Shirman, Phys. Rev. D 76, 115002 (2007) [ar Xiv:0706.2677l [hep-ph]]: 
T. G. Rizzo, JHEP 0710, 044 (2007) jarXiv:0706.3025l [hep-ph]]; H. Goldberg and P. Nath, 
larXiv:0706.3898ll [hep-ph]; L. Anchordoqui and H. Goldberg, Phys. Lett. B 659, 345 (2008) 
|arXiv:0709.0678] [hep-ph]]. R. Zwicky, larXiv:0707.0677] [hep-ph] ; larXiv:0710.4 430] [hep-ph]; 
T. Kikuchi and N. Okada, larXiv:0707.0893] [hep-ph]; larXiv:071 1.1506] [hep-ph]; R. Mohanta 
and A. K. Giri, Phys. Rev. D 76, 057701 (2007) [arXiv:0707.3308] [hep-ph]]: larXiv:0711.3516l 
[hep-ph]; C. S. Huang and X. H. Wu. larXiv:0707.1268l [hep-ph]: N. V. Krasnikov, Int. J. Mod. 
Phys. A 22, 5117 (2007) [arXiv: 0707. 1419] [hep-ph]]; H. Zhang, C. S. Li and Z. Li, Phys. 



16 



Rev. D 76, 116003 (2007) |arXiv:0707.2T32] [hep-oh]]. Y. Nakayama, Phys. Rev. D 76, 105009 
(2007) |arXiv:0707.245T] [hep-ph]]; N. G. Deshpande, X. G. He and J. Jiang, Phys. Lett. B 
656, 91 (2007) [arXiv: 0707. 2959) [hep-ph]]; N. G. Deshpande, S. D. H. Hsu and J. Jiang, 
larXiv:0708.2735] [hep-ph]: A. Delgado, J. R. Espinosa and M. Quiros, JHEP 0710, 094 (2007) 
|arXiv:0707.4309] [hep-ph]]: M. Neubert. [i;SCiv:0708.0036l [hep-ph]: S. Hannestad, G. Raffelt 
and Y. Y. Y. Wong, Phys. Rev. D 76, 121701 (2007) [arXiv: 0708. 14041 [hep-ph]]: P. K. Das, 
Phys. Rev. D 76, 123012 (2007) |arXiv:0708.28T2] [hep-ph]]: D. Maiumdar. rarXiv:0708.3485] 
[hep-ph]; A. T. Alan and N. K. Pak, larXiv:0708.3802] [hep-ph]; A. T. Alan, N. K. Pak and 
A. Senol, larXiv:0710.4239] [hep-ph]; A. T. Alan, larXiv:071 1.3272] [hep-ph]; A. Freitas and 
D. Wyler, larXiv:0708.4339] [hep-ph]; I. Gogoladze, N. Okada and Q. Shafi, Phys. Lett. B 
659, 357 (2008) [arXiv:0708.4405] [hep-ph]]; T. i. Hur, P. Ko and X. H. Wu, Phys. Rev. D 
76, 096008 (2007) [arXiv:0709.0 629l [hep-ph]]; S. Majhi, larXiv:0709.1960] [hep-ph]; J. Mc- 
Donald, larXiv:0709.2350] [hep-ph]; S. Das, S. Mohanty and K. Rao, larXiv: 0709. 2583) [hep- 
ph]; A. Shomer. [a?Xiv:0709.3555] [hep-th]: A. Kobakhidze, Phys. Rev. D 76, 097701 (2007) 
|arXiv:0709.3782] [hep-ph]]. A. B. Balantekin and K. O. Ozansoy, Phys. Rev. D 76, 095014 
(2007) [arXiv:0710.0028l [hep-ph]]; X. Liu, H. W. Ke, Q. P. Qiao, Z. T. Wei and X. Q. Li, 
larXiv:0710.2600l [hep-ph]; E. O. Iltan, larXiv:0710.2677] [hep-ph]; larXiv: 071 1.2744] [hep-ph] ; 
J. P. Lee, larXiv:0710.2797l [hep-ph]; G. F. Giudice, larXiv:0710.329l [hep-ph]. I. Lewis, 
larXiv:0710. 41471 [hep-ph]: G. L. Alberghi, A. Y. Kamenshchik, A. Tronconi, G. P. Vacca and 
G. Venturi. [a?Xiv:0710.4275l [hep-th] : D. I. Kazakov and G. S. Vartanov. larXiv:0710.4889] [hep- 
ph]; G. W. S. Hou. larXiv:0710.5424l [hep-ex]: O. Cakir and K. O. Ozansov. Ia?Xiv:0710.5773l 
[hep-ph]; I. Sahin and B. Sahin, larXiv : 0711 • 1665] [hep-ph]; T. A. Ryttov and F. Sannino, 
larXiv:0711.3745l [hep-th]; K. Huitu and S. K. Rai, .arXiv:0711.4754l [hep-ph]; S. Dutta and 
A. Goyal, larXiv:0712.0"T45] [hep-ph]; J. R. Mureika, larXiv:0712.1786i [hep-ph]; T. Han, Z. Si, 
K. M. Zurek and M. J. Strassler, larXiv:0712."204n [hep-ph]; C. Germani and A. Schelpe, 
larXiv: 071 2. 2243] [hep-th]: B. Holdom. la?Xiv:0712.2379] [hep-ph]: O. Cakir and K. O. Ozansoy, 
larXiv:0712.3"8T4] [hep-ph]; Y. F. Wu, D. X. Zhang. [a?Xrv:0712.3923 l [hep-ph]: T. Kikuchi, N. 
Okada and M. Takeuchi. larXiv:0801.00"T8] [hep-ph] : X. G. He and S. Pakvasa. ,iarXiv:0801.0189] 
[hep-ph]; E. O. Iltan. larXiv:0801.0301l [hep-ph]. 
[16] R. Mohanta and A. K. Giri, Phys. Rev. D 76, 075015 (2007) [arXiv: 0707. 1234] [hep-ph]]; 
A. Lenz, Phys. Rev. D 76, 065006 (2007) [arXiv: 0707. 1535] [hep-ph]]. 



17 



A. G. Akeroyd et al, |arXiv:hep-ex/040607T 



M. Bona et al. larXiv:0709.0451l [heo-ex]. 



M. Calvi, arXiv:hep-ex/0506046 



T. Banks and A. Zaks, Nucl. Phys. B196, 189 (1982). 

B. Grinstein, K. Intriligator and I. Z. Rothstein. larXiv:0801.1140] [hep-ph]. 

H. Georgi and Y. Kats, larXiv:0805.3953] [hep-ph]. 

J. Walder [DO Cohaboration], [axXiYi071M235J [hep-ex]. 

A. Abulencia et al. [CDF CoUaboration], Phys. Rev. Lett. 97, 242003 (2006) 



[arXiv:hep-ex/0609040| . 

G. Buchaha, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996) 



|arXiv:hep-ph/9512380l . 

T. Inami and C. S. Lim, Prog. Theor. Phys. 65, 297 (1981) [Erratum-ibid. 65, 1772 (1981)]. 
P. Bah and R. Fleischer, Eur. Phys. J. C 48, 413 (2006) [arXiv:hep-ph/0604249| . 
S. Aoki et al. [JLQCD Collaboration], Phys. Rev. Lett. 91, 212001 (2003) 
|arXiv:hep-ph/0307039] . 

A. Gray et al. [HPQCD Collaboration], Phys. Rev. Lett. 95, 212001 (2005) 
|arXiv:hep-lat/0507015| . 

L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). 



Unitarity Triangle fit (UTfit), http://www.utfit.org 



O. Buchmuller and H. Flacher, Phys. Rev. D 73, 073008 (2006) |arXiv:hep-ph/0507253 



updated value is available at http://www.slac.stanford.edu/xorg/hfag/ 



S. L. Chen, X. G. He, X. Q. Li, H. C. Tsai and Z. T. Wei. larXiv:0710.3663l [heo-ph]. 
M. Gronau, J. L. Rosner and D. London, Phys. Rev. Lett. 73, 21 (1994) 
[arXiv:hep-ph/9404282| ; M. Gronau, O. F. Hernandez, D. London and J. L. Rosner, 
Phys. Rev. D 50, 4529 (1994) |arXiv:hep-ph/9404 283 1 ; Phys. Rev. D 52, 6374 (1995) 
|arXiv:hep-ph/9504327l ; C. S. Kim, D. London and T. Yoshikawa, Phys. Rev. D 57, 4010 



(1998) [arXiv:hep-ph797083 56 1 . 
[35] We thanks H, n, Li et. al. for giving updated values of the topological parameters including 

the errors obtained from PQCD calculation. 
[36] W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 
[37] P. Bah and R. Zwicky, Phys. Lett. B625, 225 (2005). 



18 



[38] A. J. Buras, R. Fleischer, S. Recksiegel and F. Schwab, Phys. Rev. Lett. 92, 101804 
(2004) [ arXiv:hep-ph/0312259 ]; Nucl. Phys. B 697, 133 (2004) [ arXiv:hep-ph/04 02112|; Acta 
Phys. Polon. B 36, 2015 (2005) [arXiv:hep-ph/0410407] ; Eur. Phys. J. C 45, 701 (2006) 
|arXiv:hep-ph/0512032] ; C. S. Kim, S. Oh and Y. W. Yoon. larXiv:0707.2967l [hep-ph]. 



19 



